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The nonequilibrium dynamics of quantum fields is an initial- value problem, which can be described 
by Kadanoff-Baym equations. Typically, and in particular when numerical solutions are demanded, 
these Kadanoff-Baym equations are restricted to Gaussian initial states. However, physical initial 
states are non-Gaussian correlated initial states. In particular, renormalizability requires the initial 
state to feature n-point correlations that asymptotically agree with the vacuum correlations at short 
distances. In order to identify physical nonequilibrium initial states, it is therefore a precondition to 
describe the vacuum correlations of the interacting theory within the nonequilibrium framework. In 
this paper, Kadanoff-Baym equations for non-Gaussian correlated initial states describing vacuum 
and thermal equilibrium are derived from the 2PI effective action. A diagrammatic method for 
the explicit construction of vacuum and thermal initial correlations from the 2PI effective action is 
provided. We present numerical solutions of Kadanoff-Baym equations for a real scalar & 4 quantum 
field theory which take the thermal initial 4-point correlation as the leading non-Gaussian correction 
into account. We find that this minimal non-Gaussian initial condition yields an approximation to 
the complete equilibrium initial state that is quantitatively and qualitatively significantly improved 
as compared to Gaussian initial states. 
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I. INTRODUCTION 

Nonequilibrium processes within astro-particle and 
high-energy physics, like reheating after inflation, baryo- 
genesis, or relativistic heavy ion collisions, are typically 
described by classical or semi-classical equations. These 
include Boltzmann equations, hydrodynamic transport 
equations or effective equations of motion for a coherent 
scalar field expectation value [1-3]. The semi-classical 
treatment provides the possibility to relate fundamental 
parameters of the underlying theory with model predic- 
tions. Although inflation and baryogenesis occur at ex- 
tremely high energies, key observables like the baryon 
asymmetry and the spectral index are subject to experi- 
mental verification, for example by measurements of the 
cosmic microwave background radiation [4]. Therefore 
it is of great importance to assess the reliability of the 
underlying semi-classical approximations. This can be 
achieved by a comparison with a completely quantum 
field theoretical treatment. 

In recent years it has been demonstrated that the time 
evolution of relativistic scalar and fermionic quantum 
fields far from equilibrium can be described based on 
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first principles by Kadanoff-Baym equations [5-9]. These 
equations for the complete one- and two-point correla- 
tion functions can be obtained from the stationarity con- 
ditions of the 2PI effective action [10] defined on the 
Schwinger-Keldysh closed real-time contour [11-15]. The 
advantages of this approach are manifold: First, its con- 
ceptual simplicity is very attractive. The only assump- 
tion entering the derivation of Kadanoff-Baym equations 
is the truncation of the so-called 2PI functional, which 
amounts to a controlled approximation in the coupling 
constant or the inverse number of field degrees of free- 
dom for specific quantum field theories [6]. Furthermore, 
Kadanoff-Baym equations inherently incorporate typical 
quantum (e.g. off-shell) effects as well as 'classical' (e.g. 
on-shell) effects in a unified manner, and can be applied 
even to systems far from equilibrium. Accordingly, they 
are very versatile and can be employed both to assess 
the validity of conventional semi-classical approximations 
(e.g. for baryogenesis and leptogenesis) , and in situations 
where a single effective description does not exist (e.g. for 
(p)reheating by inflaton decay and subsequent thermal- 
ization) [16-20]. 

It has been shown that numerical solutions of Kada- 
noff-Baym equations not only provide a description of the 
quantum thermalization process of relativistic quantum 
fields for closed systems [5, 21], but also feature a sepa- 
ration of time-scales between kinetic and chemical equili- 
bration (prethermalization) [22]. Furthermore, they have 



been compared to semi-classical transport equations for 
bosonic and fermionic systems [7, 9, 23-26] (see also 
Refs. [27-31] for the non-relativistic case). Moreover, 
Kadanoff-Baym equations can describe the decay of a 
coherent, oscillating scalar field expectation value under 
conditions where parametric resonance occurs [16], and 
have also been investigated in curved space-time [32, 33]. 

These successes of the 2PI effective action and Kada- 
noff-Baym equations in the area of noncquilibrium quan- 
tum field theory make it worthwhile and, in view of re- 
alistic applications, necessary to answer remaining con- 
ceptual questions, like renormalization. The renormal- 
ization of the 2P1 effective action in vacuum and at finite 
temperature has been established recently [34-38]. ft 
has been shown that the vacuum counterterms are also 
sufficient at finite temperature. In order to extend this 
proof to nonequilibrium situations, it is necessary to iden- 
tify initial states that are themselves free of divergences. 
In particular, this requires that the correlation functions 
characterizing these physical initial states are rendered 
finite by the vacuum counterterms. 

Typically, Kadanoff-Baym equations are solved for 
Gaussian initial states. All connected n-point correlation 
functions with n > 2 vanish for Gaussian initial states by 
definition. However, in vacuum and at finite temperature 
the 3- and 4-point correlation functions carry overall di- 
vergences that are cancelled by corresponding vacuum 
counterterms. Thus Gaussian initial states lead to an 
unbalanced divergence at the initial time [39]. In or- 
der to overcome this shortcoming, physical initial states 
have to carry non-Gaussian initial 3- and 4-point corre- 
lation functions that differ from the vacuum correlations 
at most by a finite amount. 

In order to be able to identify physical initial states, it 
is therefore a precondition to be able to describe vacuum 
and thermal equilibrium within the standard framework 
of noncquilibrium quantum field theory, i.e. on the closed 
real-time path with finite initial time ti n u = 0. Apart 
from the question of renormalization it is also a matter 
of principle that vacuum and thermal equilibrium should 
be accessible within nonequilibrium field theory as special 
cases by choosing the initial state appropriately. 

In this work, Kadanoff-Baym equations for non-Gaus- 
sian correlated initial states describing vacuum and ther- 
mal equilibrium are derived from the 2PI effective action 
formulated on the closed real-time path with finite ini- 
tial time. For that purpose, we propose a diagrammatic 
method for the explicit construction of vacuum and ther- 
mal initial correlations that is applicable to nonpertur- 
bative 2PI approximations. 

There exist several techniques to describe non-Gaus- 
sian correlated initial states. These can be divided into 
two categories: Either, the correlations are generated 
by modifying the closed real-time path C, or the initial 
state is explicitly described by its density matrix p. We 
shall refer to these as implicit and explicit techniques, re- 
spectively. The implicit techniques include the so-called 
imaginary-time stepping [15, 27, 28], where an imaginary 



branch is added to the contour C, similar to the descrip- 
tion of thermal field theory [40-42]. Another possibility 
is to extend the closed real-time contour C over the com- 
plete real axis, such that it runs from — oo to +oo, and 
back to — oo. The correlated initial state is then gener- 
ated by including an external two-point source K(x,y), 
that is switched off at the 'initial' time [39]. 

For the explicit technique, the density matrix p of the 
initial state is parameterized by initial n-point correla- 
tion functions a n (x\, . . . ,x n ) [43, 44]. These appear in 
the form of non-local effective n-point vertices in the 2PI 
effective action. When deriving Kadanoff-Baym equa- 
tions, these inherit the contributions from non-Gaussian 
initial correlations. 

The advantage of the implicit techniques is that the 
equilibrium limit can be approached without any addi- 
tional work. However, due to the implicit preparation of 
the initial state, the freedom and the control in choos- 
ing the initial state is restricted. The advantage of the 
explicit technique is that the resulting Kadanoff-Baym 
equations are very similar to the Gaussian case. Further- 
more, the explicit approach provides a maximal degree 
of freedom for specifying the initial state. 

The renormalization of Kadanoff-Baym equations has 
recently been discussed based on the implicit technique 
involving an external two-point source [39] . In this work, 
we use the explicit technique. Thus, the methods devel- 
oped in this paper provide a complementary framework 
for addressing the issue of renormalization. 

The main purpose of this work is to provide tech- 
niques for calculating the non-Gaussian initial correla- 
tions cn„ (x\, . . . , x n ) for a thermal initial state p t h within 
the 2PI-Schwinger-Kcldysh formalism. As explained 
above, this is a prerequisite for studying the renormal- 
ization of Kadanoff-Baym equations based on the explicit 
approach. In Ref. [44], a perturbative expansion of the 
ct 1 ^ was derived. Unfortunately, this expansion is not 
suitable for the non-perturbative 2PI formalism. The 
main idea followed in the present work is to determine 
the cn„ by matching the Kadanoff-Baym equations for 
a thermal initial state on the one hand with the evo- 
lution equations obtained from the 2PI effective action 
formulated on the well-known thermal time contour [40- 
42, 45, 46] on the other hand. For this matching pro- 
cedure, we employ the thermal time contour obtained 
by concatenating the closed real-time contour with the 
imaginary time contour. We stress that it is important 
to keep a finite 'initial' time. In this case, both the hori- 
zontal and the vertical branches contribute [45-47]. The 
contributions from the vertical branch can then be identi- 
fied with corresponding contributions from the a„ within 
the equivalent Schwinger-Keldysh formalism. We note 
that, when considering an exact thermal initial state, the 
propagator obtained from the Kadanoff-Baym equations 
is time-translation invariant within its domain of defini- 
tion, in accordance with Refs. [45-47] . 

This work is organized as follows: In section II, we 
derive Kadanoff-Baym equations for non-Gaussian ini- 
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FIG. 1: Closed real-time path C. This time path was invented 
by Schwinger [11] and applied to non- equilibrium problems by 
Keldysh [14] (see also Refs. [12, 13]). In order to avoid the 
doubling of the degrees of freedom, we use the form presented 
in Ref. [15]. 

tial states from the corresponding 2PI effective action 
using the explicit technique. In section III, we pro- 
vide techniques for calculating the non-Gaussian initial 
correlations a^(x\, . . . , x n ) within perturbation theory. 
These techniques are generalized to the nonperturba- 
tive 2PI case in section IV. In section V, Kadanoff- 
Baym equations for a thermal initial state are derived 
by combining the results from sections II and IV. In 
section VI, we compare numerical solutions of Kadanoff- 
Baym equations for two noncquilibrium initial states that 
are approximations to the thermal initial state. One of 
them is Gaussian, while the other also includes the lead- 
ing non-Gaussian initial correlation. The appendices A 
and B contain additional material helpful for sections III 
and IV, respectively. 

II. EFFECTIVE ACTION AND 
KADANOFF-BAYM EQUATIONS 

A. Gaussian Initial States 

In this subsection, we start from the classical action for 
a real scalar quantum field with a quartic self interaction 

S[<t>] = J d*x Q(90) 2 - ± m v - (1) 

and review the basic elements of the derivation of the 
2PI effective action and the Kadanoff-Baym equations 
for the case of a Gaussian initial state. In the following 
subsections, we can then easily expose the differences, 
which arise for a non-Gaussian initial state. 

The Schwinger-Keldysh propagator is defined by 

G(x, y) = {T C $(aO*(y)> - (*(*)> <<%)> , (2) 

where Tc denotes the time-ordering operator along the 
closed real-time path C shown in figure 1 [11, 14, 15]. 
The Schwinger-Keldysh propagator can be obtained by 
functional differentiation from the generating functional 
for correlation functions formulated on the closed real- 
time path. The generating functional in the presence of 



a local external source J (a;) and a bilocal external source 
K(x,y) is given by [44] 

Z P [J,K] = J (x) (<p+\p\<p-) 

x exp ^iS [ip] + iJtp + ^ipKifij , (3) 

where a matrix-vector notation has been used for the 
space-time integrals in the exponential function 1 , and 
\ip±) are the quantum states corresponding to the field 
configurations ip±(x) = ip(0±,x). The information 
about the initial state enters via the matrix element of 
the density matrix p, which is known only at the initial 
time t = t lnit = 0. 

A Gaussian initial state is an initial state for which 
all connected n-point correlation functions with n > 3 
vanish at the initial time. The density matrix clement 
for a Gaussian initial state can be parameterized by 

(tp + \p\tp-)=exp(ia + ioti<p+^<pa2<pj ■ (4) 

Therefore, in the Gaussian case, the contribution of the 
density matrix to the generating functional (3) can for- 
mally be absorbed into the external sources, J + a\ — > J 
and K + a 2 ~ > K. As is, for example, explained in 
Ref. [48], this means that the Gaussian initial density 
matrix does not appear explicitly, but rather enters the 
dynamics via defining the initial conditions of all in- 
dependent one- and two-point functions. For vanish- 
ing field expectation value, these are given by G(x,y), 
(d x o + d y o)G(x,y) and d x ad y oG(x,y), all evaluated at 
x° = y° = [9, 48]. 

The 2PI effective action T[d>, G] is the double Legen- 
dre transform of the generating functional (3) with re- 
spect to the external sources. For a Gaussian initial 
state, the generating functional has the same structure 
as the generating functional in vacuum, except that all 
time- integrations are performed over the closed real-time 
path. Consequently, for a Gaussian initial state the 2PI 
effective action can be parameterized in the form [10] 

T[4>,G\ = SM + ^Trlog^G- 1 ] 

+ ^Tr c [g - 1 G]+r 2 [^,G] , (5) 

where Go 1 is the inverse classical Schwinger-Keldysh 
propagator and iT2[<p,G] is the sum of all 2PI Feyn- 
man diagrams without any external legs, where internal 
lines represent the complete Schwinger-Keldysh propa- 
gator G(x,y). The vertices of the diagrams contained 



Throughout this work, the compact notation of Ref. [15] is used 
for the contour integrals over the closed real-time path, for exam- 
ple Jf = J c d 4 x J(x)ip(x) = f d 4 x \Jj r {x)ip+(x) — J-(x)ip-(x)]. 
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in iT2 [4>, G] are given by the third and fourth functional 
derivatives of the classical action S[4>] [10]. Eventually, 
the Kadanoff-Baym equations 

v° 

(a x +M 2 (x))G F (x,y) = J d i zU F (x,z)G p (z,y) 

o 

- J d 4 zU p (x,z)G F (z,y) , (6) 
o 

v° 

(n x + M 2 (x))G p (x,y) = f d 4 zTl p (x,z)G p (z,y) , 



follow from the stationarity condition of the 2PI effective 
action. Here, we use the notation of Ref. [9]. 



B. Non-Gaussian Initial States 

In the remainder of this section, Kadanoff-Baym equa- 
tions are derived that can describe systems characterized 
by a general non-Gaussian initial state [43, 44, 49]. For 
that purpose, we extend the derivation of the previous 
subsection using a generalization of the Gaussian density 
matrix (4). In general, the matrix element of the density 
matrix p is an arbitrary functional of the field configura- 
tions <p+{x) and <p-(x), which can be written as [44] 



(P+\p\P-) =exp(iF[<^]) . 



(7) 



While for a Gaussian initial state F[<p] is a quadratic 
functional of the field, for a general non-Gaussian initial 
state it may be Taylor expanded in the form [44] 



F[p] 



E 

n=0 



— / d i x 1 . . . d A x n a n (xi , . . . , x n ) 
n\ J 



x tp(xi) 



p(x n ) 



(8) 



By definition F[<p] depends only on the field configuration 
evaluated at the boundaries of the time contour. Con- 
sequently, the kernels a n {x\, . . . , x n ) are non-zero only 
if all their time arguments lie on the boundaries of the 
time contour. With the notation 5+(t) = Sc(t — 0+) and 
5-(t) = 5c(t — 0_), they can be written in the form 

&n {,X\ i • • • > ^n) — &n ' ' ™ ) • • • ) *^n) 

x 6 €1 (*») • . . . • <J e „ (x° n ) , (9) 

where summation over ej G {+, ~ } is implied. In this 
way, the explicit dependence of the functional F[ip] on the 
field configurations ip+(x) and ip~(x) may be recovered, 

F[ip] =a + J d 3 x a{(x)<p e {x) + ... . 

The set of all kernels a n with n > encodes the complete 
information about the density matrix characterizing the 



initial state. Not all the kernels are independent. The 
Hermiticity of the density matrix, p = p\ implies that 



,asn) = (*Q4- ei )--(- e »)(xi,...,x n ))* . 



If the initial state is invariant under some symmetries, 
there are further constraints. For example, for an ini- 
tial state which is invariant under the Z 2 -symmetry 
<f> — > — $, all kernels a n (x\, . . . , x n ) with odd n van- 
ish. If the initial state is homogeneous in space, the ini- 
tial correlations a n (xi, . . . , x n ) are invariant under space- 
translations Xi — > Xi + a of all arguments for any real 
three- vector a, and can conveniently be expressed in spa- 
tial momentum space, 



a 



ei,...,e r 



{xi , • ■ • , X n ) 



d 3 k x 



^ exp t 



(2tt) 3 "J (2tt) 3 
x (2n) 3 5 3 (k 1 + --- + k n )a%'"" e " (fc 1; . . . , . 
Summarizing, the generating functional reads 

Z P [J,K] = j V P exp(iS[<p] 

+ iJ(p+^tpK(p + iF 3 [tp]) , (10) 

where the kernels ot\ and a 2 have been absorbed into 
the measure T>ip and into the sources J and K, respec- 
tively. The functional F 3 [ip] = F 3 [ip; a 3 , a^, . . . ] contains 
the contributions of third, fourth and higher orders of the 
Taylor expansion (8) and vanishes for a Gaussian initial 
state. 



C. 2PI Effective Action for Non-Gaussian Initial 
States 

According to eq. (10), the 2PI effective action in the 
presence of non-Gaussian correlations is obtained from 
the standard parameterization [10] of the 2PI effective 
action applied to a theory described by the modified clas- 
sical action S[<f>] = S[<j>] + F 3 [<j>], 



T[0,G] = ^R-TrlogcfG- 1 ] 



+ 2 Trc 



g 1 c\ + f 2 [<j>,G] 

= T G [(j>,G]+T nG [ct>,G] , 



(11) 



where ^{x.y) = 6 2 S[cp}/5(j)(x)8(f>( y ). The Gaus- 
sian part r<3 [4>, G] coincides with the right-hand side of 
eq. (5), and the non-Gaussian part is given by 



T nG [d>,G]=F 3 [ ( f>} + -Tt c 



5 2 F 3 



G 



T2,nG[0, G] 
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FIG. 2: Non-local effective vertices ia n (x\, . . . , x n ) connecting 
n lines for n = 3,4,5,6 encoding the non-Gaussian three-, 
four-, five-, and six-point correlations of the initial state. 



FIG. 3: Diagrams contributing to the three-loop truncation of 
the 2PI effective action in the symmetric phase (setting-sun 
approximation) in the presence of an effective non-local four- 
point vertex. 



where 



r 2 , nG [^,G] = f 2 [^G]-r 2 [0,G] . 



The modified 2PI functional iT 2 [cf), G] is equal to the sum 
of all 2PI Feynman diagrams without any external legs, 
where internal lines represent the complete propagator 
G(x, y) and where vertices are given by the functional 
derivatives of the modified classical action S[4>] — S[<j>] + 
F 3 [(f>]. The contribution from the classical action S[<f>] 
leads to the classical local three- and four-point vertices. 
Additionally, the contribution from the functional ^[0] 
leads to effective non-local vertices, which contain the 
non-Gaussian initial n-point correlations with n > 3 (see 
figure 2), 



8<j)(xi) . . . 84>(x n ) 



— iOi n {x\ , . . . , Xn) • 



(12) 



These effective n-point vertices are only supported at 
the initial time, and can be parameterized analogously 
to eq. (9). For a Z 2 -symmetric initial state, the 
field expectation value vanishes, <fi(x) = 0, such that 
a n (xi, . . . , x n ) — a n (xi, . . . , x n ). The contribution of 
these effective non-local vertices is most important close 
to the initial time. For example, a non-zero four-point 
source a.i(xi,. .. ,x±) leads to a non- vanishing value of 
the connected four-point correlation function at the ini- 
tial time, which is impossible for a Gaussian initial state. 

Note that those 2PI diagrams that contain exclu- 
sively the classical vertices contribute to the functional 
ir 2 [0, G] by definition. Therefore, the diagrams con- 
tributing to the non-Gaussian part iT 2tn c[(j) 1 G] contain 
at least one effective vertex from cq. (12). 

In section VI we study the numerical solution of Ka- 
danoff-Baym equations for a Z 2 -symmetric non-Gaussian 
initial state with a non-zero initial 4-point correlation. In 
this case the 2PI functional T 2 reads in "naive" 2 thrce- 



2 This means that non-local effective vertices do not affect the 
counting of loops. 



loop approximation (see figure 3) 
1 



«r 2 [G] 



2^1234 



+ 



48/ 



d x\. 



X Gi5G 2 6G37G48 



- «A<5i 2 <5 2 3<534 + iai234 
zA<5i 2 5 2 3<534 + iai 2 34 

- iX5 56 S 57 8 5S + ««5678 



Gi 2 G34 

(13) 



where G i2 = G(xi,x 2 ) and ai 23 4 = a^xi, x 2 , X3, X4). 
Note that the contribution to the mixed basketball dia- 
gram in the second and third line with one classical and 
one effective vertex appears twice, which accounts for the 
symmetry factor 1/24. 



D. Self-Energy for Non-Gaussian Initial States 

The equation of motion for the complete propagator 
obtained from cq. (11) reads 

G^{x,y) = g^ 1 (x 7 y) - U(x,y) - ia 2 (x,y) , (14) 

where <5 2 = a 2 + S 2 F^[(j)]/S(f)8(j) and the complete self- 
energy is given by 



U(x,y) = 



2iST 2 [(j),G} 

5G(y,x) 
2i5T 2 [4>,G] , 2i5T 2 , nG [ 



,G] 



5G(y,x) SG(y,x) 
= U G {x,y)+n nG (x,y), (15) 

where II contains the contributions to the self energy, 
which are also present for a Gaussian initial state, and 
the non-Gaussian part U nG contains diagrams with at 
least one non-local effective vertex. They can be further 
decomposed as 

n G (a;,2/) = -iUi oc (x)5 c (x-y) + U G on _ loc (x,y), 
U nG (x,y) = zU: G rface (x,y) + U n n G n _ loc (x,y). (16) 

The non-Gaussian non-local part n™^_ (oc (x, y) con- 
tains diagrams where both external lines are attached 
to a standard vertex. The non-Gaussian "surface" part 
^™urface( x > v) contains diagrams where at least one ex- 
ternal line is attached to a non-local effective vertex. 
Thus, the surface part is supported only at the initial 
time surface where x° = or y° = 0. In general, such 
contributions can arise in the following ways: 
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FIG. 4: Diagrams contributing to the self-energy H(x,y) in 
setting-sun approximation in the presence of an effective non- 
local four-point vertex. From left to right, the diagrams con- 



tribute to II/, 



n 



and the last two both 



contribute to H aa - 



1. From diagrams where both external lines are con- 
nected to an effective non-local vertex as given in 
eq. (12). They arc supported at x° = y° = 0. 

2. From diagrams where one of the two external lines 
is connected to an effective non-local vertex, while 
the other one is connected to a classical local vertex. 
They are supported at x° — 0, y° > or vice- versa. 

3. Via the contribution ia 2 (x,y) of the initial two- 
point source which is supported &t x° = y a = 0. 
This is the only Gaussian surface-contribution. 

Accordingly, the contributions to the self-energy which 
are supported at the initial time surface can be further 
decomposed as 



^-surface {x, y) — ^-surface 



= n n 



(x,y) + a 2 {x,y) 
(x, y) + n AQ (x,y)+ n aX (x, y) , 



where 



U aa (x,y) 
K\a{x,y) 
R a \(x,y) 



IL\ a (x°,x,y)6 e (y ), (17) 
6 e (x°)W aX (x,y°,y) = U Xa (y,x). 



Tl aa contains all contributions of type (1.) and (3.). Di- 
agrams of type (2.) contribute to U Xa or IT qA depending 
which external line is attached to the effective non-local 
vertex and which to the classical local vertex. For all di- 
agrams contributing to H Xa the left line is connected to 
the classical four- or three-point vertex. The non-local 
part of the self-energy can be split into statistical and 
spectral components, similarly to the Gaussian case, 



^-non — loc{x, y) — ^-non—lni ' 



(x,y)+TV^ n _ l0C (x,y) (18) 



= U F (x,y) - -sign c (x° - y°)IL p (x,y) . 



The local part of the self-energy is identical to the Gaus- 
sian case and is included in the effective time-dependent 
mass term M 2 (x) = m 2 + Ui oc (x). 

For the setting-sun approximation from eq. (13), the 
self-energy is given by (see figure 4) 



IT 



non—loc 



n 



nG 



non — loc x ■ ~ , '1 /" 

m aa (x,y) = ia 2 (x,y) + - d 4 x 34 ia xy34: G 3 4: 

+ 6 j rf4x2 "' 7 * aa; 234G25G 36 G , 47 ia 567 y , 

til\a(x,y) = Jd 4 x 123 G xl G x2 G x3 ia 123y , 



{x,y) = 

{x,y) 



G{x,yf 



iU aX (x,y) = Jd A x 23i ia x23i G 2y G 3y G iy . 



E. Kadanoff-Baym Equations for Non-Gaussian 
Initial States 



Convoluting eq. (14) with the complete propagator 
yields 



(a x + M 2 {x)) G(x, y) = -i5 c {x - y) 



(19) 



■; J d 4 z [n 



non — loc 



(x, z) + iU Xa (x, z)]G(z,y) . 



The second line follows from the parameterization (16) 
of the self-energy, and assuming x° > and y° > 
0. Using eqs. (17,18) and transforming to spatial mo- 
mentum space yields the Kadanoff-Baym equations for 
Gf(x° ,y° ,fe) and G p (x°,y°,k) for spatially homoge- 
neous non-Gaussian initial states, 



y x- 
(d 2 x0 +k 2 + M 2 (x )) G F (x°,y°, k) = Jdz° n F (x°, z°, k)G p (z°, y a ,k) - Jdz° U p (x°, z°, k)G F (z°, y°, k) 



+ n AQ , F (x°, k)G F {0, y°, k) + -n Aa!( ,(x°, k)G p (0, y°,k) 



(20) 



and 



y 

(d 2 x o+k 2 + M 2 (x ))G p (x°,y°,k) = Jdz n p (x ,z ,k)G p (z 7 y°,k). 







c 












— : ^- 


0_ 




tinit 






X 


-i0 







FIG. 5: The thermal time contour C+X is obtained by concate- 
nating the closed real-time contour C and the imaginary time 
contour X running from t = to t = —if) [40-42, 45, 46]- 
We employ this time path in order to infer the initial correla- 
tion functions a„(xi, . . . , x„) required for describing thermal 
equilibrium on the Schwinger-Keldysh closed real-time path C 
with finite initial time shown in figure 1. 



where 

= n+ a (x\k) + xr Xa (x\k), 
iWx°,fc) = 2i(n+ a (x°,k)-n^ a (x°,k)) . 

Thus, for a non-Gaussian initial state, the right-hand 
side of the Kadanoff-Baym equation for the statistical 
propagator is modified. In addition to the memory in- 
tegrals there are now new contributions originating from 
the non-Gaussian initial correlations. Unlike the memory 
integrals, these new contributions do not have to vanish 
in the limit x°,y° — > 0. This is due to the fact that the 
higher non-Gaussian correlations of the initial state can 
lead to a non- vanishing value of the connected four- and 
three-point correlation functions at the initial time. 



III. THERMAL INITIAL CORRELATIONS: 
PERTURBATION THEORY 

In order to derive Kadanoff-Baym equations that are 
capable of describing thermal equilibrium, the thermal 
density matrix 

Pth = ^ exp {-0H) , 

has to be represented by a Taylor expansion in terms 
of thermal correlation functions a„ (ari, . . . , x n ) as in 
eqs. (7,8). These thermal correlation functions do then 



enter the Kadanoff-Baym equations in the form of non- 
local effective vertices, as described in the previous sec- 
tion. 

The thermal correlations functions a^(x\, . . . , x n ) can 
be calculated order-by-order in the coupling constant 
within usual perturbation theory (see appendix A). How- 
ever, in the context of Kadanoff-Baym equations, it is 
necessary to use approximations of the thermal correla- 
tion functions that are compatible with the underlying 
truncation of the 2PI effective action. It is a major pur- 
pose of this paper to provide computational techniques 
for identifying suitable approximations of the thermal 
correlation functions a%(xi, . . . , x n ). 

For simplicity, in this section, we first present the com- 
putational techniques within perturbation theory. In the 
following section, these techniques are then generalized 
to the nonperturbative 2PI formalism. 

The main idea is to determine the functions 
a^(xi, . . . ,x n ) by matching the description of thermal 
equilibrium based on the closed real-time path C in the 
presence of effective vertices a* ("C + a") on the one 
hand with the well-known equivalent description based 
on the thermal time path (" C+X") [40-42, 45, 46] shown 
in figure 5 on the other hand. 

The generating functional for these two descriptions 
is obtained by inserting the respective representations of 
the thermal density matrix, 



' <p(-i(),x)=<p + (x) 

j T>ip exp 

<p(0,x)=ip- (as) 



(f+ \(Hh\f-} 



(i J d 4 xC(x) 



for "C+X" , 



exp 



Eth 
"12- 



* > "l2...nVlV2 ■■■<Pn 
n=0 



for "C + a' 



(21) 

into cq. (3). The argument of the exponential in the lower 
expression is a short-hand notation for eqs. (7, 8). In the 
following, we show how perturbative Feynman diagrams 
formulated within the well-known C+I formalism can 
equivalently be represented within the C + a formalism. 



A. Thermal time contour C+X 

In this subsection, we briefly review the well-known 
description of thermal equilibrium based on the thermal 
time contour C+X in order to establish the notations 
required later on. We stress that, for our purpose, we 
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have to keep the initial time finite. For this case, the 
formulation of thermal field theory has been discussed in 
Refs. [45-47]. The free thermal propagator defined on 
the thermal time contour C+T is 

iG o*th( x > V) = (- u * - ™ 2 ) 6 c+x(x - y) , 

for x° ,y° e C+T. It may be decomposed into the free 
thermal statistical propagator Go.F(x,y) and the free 
thermal spectral function Go, p (x,y), 



We note that So(x,y)\ x o— y o is time-independent, as ex- 
pected in thermal equilibrium. Nevertheless, the four 
contributions shown in the last line may individually de- 
pend on time. However, this time-dependence cancels 
out in their sum, as has been, for example, discussed in 
Rcf. [46]. 



Closed real-time contour with thermal initial 
correlations C + a 



G .th(x,y) = G 0:F (x,y) - -sign c+x (a; - y°)G , p (x, y) . 
The explicit solution of the free equation of motion is 

G 0tF (x°,y l \k) = cos (uj k (x ~ y°)) , 

G . p (x°,y a ,k) = — sm(u k (x° -y )) , (22) 

for x°,y° E C+T. Here, riBs(wfc) is the Bose-Einstein 
distribution function, 



riBE(^k) 



1 



= \J m? + k 2 . 



Each of the two time arguments of the propagator can ei- 
ther be real or imaginary, which yields four combinations 



cc 

0,th> u 0.tti "O.t/l! 



G 



a 



Gq t h , Gq ~i h . These appear in perturbative 



Feynman diagrams which are constructed with the free 
propagator Go,th and the classical vertices. In position 
space, each internal vertex of a Feynman diagram is in- 
tegrated over the thermal time contour C+T. In order to 
disentangle the contributions from the real and the imag- 
inary branch of the time contour, the following Feynman 
rules arc defined, 



(x,y) 



r<cc 



G^ h (x,y) = o-- 



(x,y) 



r cx 

u 0,th 



G%%{x,y) 



Filled circles denote a real time, and empty circles de- 
note an imaginary time. As an example, the perturba- 
tive setting-sun diagram is considered with propagators 
attached to both external lines, and evaluated for real 
external times x°,y° € C. Both internal vertices are in- 
tegrated over the two branches C and T, respectively. 
Using the Feynman rules above, the resulting four con- 
tributions can be depicted as 

S (x,y) = -—<£—-©-—. 

( f [' 

= — g — d 4 u d 4 v G 0>t h{x,u)Gl th {u,v)G , th {v,y) 

c+x c+x 



Within the C + a formalism, all internal vertices of 
Feynman diagrams are just integrated over the closed 
real-time path C. However, the diagrams may contain 
non-local effective n-point vertices. These represent the 
n-point correlations a„ (a;i, . . . ,x n ) of the thermal initial 
state. In the following, we show how to determine the 
by a matching procedure employing the equivalent C+T 
formalism. 

Let us consider a Feynman diagram within the C+T 
formalism, like for example the perturbative setting-sun 
diagram So(x,y). We assume that all time arguments 
corresponding to the external lines are real. It turns out 
that a single diagram within the C+T formalism is rep- 
resented by a set of diagrams within the C+a formalism. 
Some of these will contain non-local effective vertices. 
Since we are working in the framework of perturbation 
theory, we have to insert approximations to the exact ef- 
fective vertices. For this purpose, we have to determine 
(i) the topologies of the required diagrams and (ii) the 
proper approximations for the effective vertices. 

In order to do so, we first consider the free thermal 
propagator evaluated with one imaginary and one real 
time. Using eq. (22) together with elementary trigono- 
metric addition theorems, it can be written as 

C IC ( ir 7/° k\ - G Mfr(~^'°' fc ) r CC / n j.) 

G 0tth { ir, y ,k) GoM0M G oM°>y > k ) 

+ id T G^ h (-ir,0,k)G^ c p (0,y°,k) . 

Next, the unequal-time statistical propagator and the 
spectral function are rewritten as 

G c o c F (0,y°,k) = Jdz S s (z°)G c c th (z°,y,k), 
c 

G c o c p (0,y°,k) = -2 t Jdz S a (z )G^ th (z ,y°,k), 



where 



S,{z°) = \{5e(z°-0 + ) + Sc(z°-0-)) , 
S a (z a ) = hsc(z -0 + )-6c(z°-0-)) . 



Combining these equations, a helpful expression for the 
free propagator evaluated with one imaginary and one 
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real time is obtained, 

G% e th (-iT,y°,k) = JdtA (-iT,t,k)G^ c th (t,y°,k), 
c 

o- ■• = o |~- -~ • . (23) 

Here, the free connection Ao(— it, z , k) is given by 
A (-i T ,z°,k) - AgHT,fc)<5 s (z°) + ASHT,fc)<5 a (z°) 
I- , (24) 

where 

<^o,th(0, U,fc) 
AS(-ir.fc) = 2d T G^ h (-ir,0,k) . (25) 

Analogously, the free propagator evaluated with one real 
and one imaginary time can be written as 

G% x th (y°,-iT,k) = JdtG% c th (x ,t,k)AT(t,-iT,k), 



(26) 



where A^(z°,-ir,k) = A (-iT,z°,k) = — ( . The 

connections A and Aj are attached to an imaginary 
and a real vertex on the left and right sides, respectively. 
Their Fourier transform into position space is 



A (v,z) 



d 3 k 
(2nf 



e ifc (»-*) A (v Q ,z°,k) , (27) 



for v° £ X and z° e C, as well as Aq(z,v) — Aq(v,z). 

In general, for any thermal diagram on C+X with V 
internal vertices, there are 2 V possibilities to combine 
the integration over C or X at each vertex. For each of 
these 2 V contributions, all lines connecting a real and an 
imaginary vertex are replaced using relations (23) and 
(26). Thereby the parts containing X- integrations are 
encapsulated into non-local effective vertices. Thus, any 
thermal diagram on C+X can equivalently be represented 
by 2 V diagrams on C, which contain classical vertices as 
well as non-local effective vertices. 

For example, the setting-sun diagram with one real and 
one imaginary vertex can be rewritten as 



-— f— -<-- 



— -f- 



(28) 

According to the symbolic notation employed here, the 
subdiagram containing the imaginary vertex, marked by 
the box, can be encapsulated into an effective non-local 
4-point vertex. Its structure is determined by the con- 
nections A and A J. This can be seen by rewriting the 



above diagrams in terms of the corresponding formal ex- 
pressions (only the first and last one are given here), 



HA) 2 
6 



d 4 u J d 4 v G 0<t h{x,u)Go tth (u,v)G ,th{v,y) 
x 

d 4 u d 4 z 1234 G ,th (x, u)G Qtth (u, zi)Go,th (% z 2 ) 



C X 

-iX 

c c 

x G .th(u,z 3 ) 



G ,th{zi,y) ■ 



In the last line, the thermal effective 4-point vertex has 
been introduced, 



a i .o L {zi,z 2 ,Z3,,z i ) 



th 
OL 



—i\ J d 4 v Ao(v, zi)Aq(v, zi) 
x 

x A {v,z 3 )A (v,z 4 ) , 

X X 





X' X 



Since the connection Aq(v, Zi) is supported only at the 
initial time zf — 0±, the effective 4-point vertex vanishes 
as soon as one of the four real times zf, . . . , z\ lies be- 
yond the initial time. Thus, the effective 4-point vertex 
has precisely the structure of a non-local effective vertex 
describing an initial correlation. Furthermore, the above 
4-point vertex constitutes the leading order contribution 
to the loop expansion of the thermal initial 4-point cor- 
relation function (see appendix A). 

Diagrams with internal lines connecting two imaginary 
vertices contain the propagator G xx h {— it, — it 1 , k). In 
order to identify the correct effective vertices in this case, 
the following relation is employed 

Go%{-iT,-iT',k) 

= D (-iT,-iT',k) + J dw° J dz° A (-iT,w°,k) 
c c 

xG c %{w\z\k)A T {z\~iT',k) 
= D Q (-iT, -it', k) 

+ A*Ht, k)G°%{0, 0, k)A s (-i T ', k) , (29) 



o o = o 



o + o ] f 



It can be verified by explicit calculation from eqs. (Al, 
22, 24). Here the propagator Dq(— it, — it', k), which is 
defined in eq. (Al), is represented by the dotted line. 
It connects two imaginary times and furnishes the per- 
turbative expansion of the thermal initial correlations as 
discussed in appendix A. Applying the upper relation, 
the setting-sun diagram with two imaginary vertices can 
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be rewritten as 



"C+X" and "C + a" can, in summary, be written as 



So(x,y) = •—-4-—-*-- 



^—4-.—^—. = ...|. ..4.....^.. ..( 

A"/, 

— • — |-» 6—| — • -(- «--|-'-p b— I — » 



- 9— -<» 



■ + ■ 



o o • - 



•6 o • 



A"/, A"/, 

y«v y-v" 
« '<""]--• + .--r*]--|* 



In the first step, the propagators connecting real and 
imaginary vertices were replaced using relation (23). 
This already yields an effective non-local two-vertex, as 
indicated in the third diagram in the first line. In order 
to check that this effective non-local two- vertex is indeed 
composed from the thermal initial correlations, the three 
propagators connecting the two imaginary vertices are 
replaced using relation (29). Accordingly, the diagram 
splits into eight terms. These can be combined to the four 
incquivalcnt contributions shown in the second line 3 . Fi- 
nally, the parts containing imaginary vertices and dotted 
lines can be identified with the corresponding contribu- 
tions to the perturbative expansion of the thermal initial 
correlations discussed in appendix A. This is represented 
graphically by encapsulating the subdiagrams inside the 
boxes. In the third line, the first diagram thus contains 
a thermal effective two-point vertex, which itself appears 
at two-loop order in the perturbative expansion of the 
thermal initial correlations. Similarly, the thermal effec- 
tive four- and six-point vertices contained in the second 
and third diagram, respectively, appear at one- and zero- 
loop order in the perturbative expansion of the thermal 
initial correlations. The two effective four-point vertices 
contained in the fourth diagram are identical to the one 
already encountered in cq. (28). 

Thus, using the representation (23) of the free propaga- 
tor connecting a real and an imaginary time and eq. (29), 
any perturbative thermal Feynman diagram formulated 
on the thermal time contour C+X can be related with a 
set of perturbative Feynman diagrams formulated on the 
closed real-time contour C. Furthermore, the required 
approximation to the complete thermal initial correla- 
tions a„ can be explicitly constructed with the help of 
the formalism introduced here. For example, for the per- 
turbative setting-sun diagram, the equivalence between 



3 Note that the symmetry factors are taken into account properly. 
For example, the symmetry factor of the second diagram in the 
second line is one third times the symmetry factor of the original 
diagram in the first line. Since there are three possibilities to 
obtain this diagram from the first one, it is obtained with the 
correct pref actor. 



+ • 



As has been noted before, So(x,y)\ x o =y o is time- 
independent. Nevertheless, the individual contributions 
shown above may depend on time, similarly as for the C+ 
X-formalism. Since Sq(x, y)\ x o =y o is time- independent, it 
is clear that the time-dependence of the individual contri- 
butions has to cancel when summing them up. Thus the 
time-translation invariance of thermal equilibrium within 
the C+a-formalism is manifestly inherited from the C+X- 
formalism [45-47] by the matching procedure described 
here. Since the same argument applies for the 2PI case 
discussed below, we will not repeat it there. 



IV. THERMAL INITIAL CORRELATIONS: 2PI 

In this section, the perturbative techniques introduced 
in the previous section are generalized to the nonpertur- 
bative 2PI case. This is required since Kadanoff-Baym 
equations are based on the 2PI formalism. While the 
Feynman diagrams shown in the previous section con- 
tain the free propagator, the diagrams treated here con- 
tain the complete propagator, which itself is the solution 
of a self-consistent Schwinger-Dyson equation. 



A. Thermal time contour C+X 

The complete thermal propagator defined on the ther- 
mal time path C+X satisfies the self-consistent Schwinger- 
Dyson equation derived from the 2PI effective action in 
thermal equilibrium, 

G- h \x,y) = i{U x + m 2 )S c+x (x - y) - U th (x, y) , (30) 

where x°,y° g C+X and Yl t h{x,y) is the thermal self- 
energy. For example, in setting-sun approximation it 
reads 

— iX (— i\) 2 

Hth{x,y) = -z-G t h(x,x)dc+x(x-y)-\ — G th (x,yf . 

2 

The complete propagator furnishes the expansion of the 
2PI effective action in terms of 2PI Feynman diagrams. 
Similar to the perturbative case, the following Feynman 
rules are defined, 



G%(x,y) = 



Gf*(x,y) = o- 



* . G%{x,y) = •- 
o , Gff(x,y) = o- 



-ixffx = x, -^jfx = X' = X 
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Accordingly, diagrams containing the complete propaga- 
tor can be decomposed in analogy to the perturbative 
case. For example, 



S(x,y) = « 
HA) 2 



J d A u J d A v G t h(x,u)G 3 th (u,v)G th (v,y) 



c+i c+i 



B. Closed real-time contour with thermal initial 
correlations C + a 

As for the perturbative case, we will now use the equiv- 
alence between C+I and C + a in order to infer the 
proper thermal initial correlations a^(x\, . . . , x n ) within 
the 2PI framework. 

In order to disentangle the contributions from the real 
and the imaginary branch of the thermal time contour, 
a generalization of eq. (23) to the 2PI case is required. 
Due to the nonperturbative nature of the 2PI formalism, 
this generalization is non-trivial and requires a somewhat 
lengthy calculation that can be found in appendix B. The 
most important result is that the complete propagator 
connecting imaginary and real times can be decomposed 
into a convolution of a complete connection A(—ir, z°, k) 
and the complete real-real propagator, 

Gl c (-ir,y°,k) = JdtA(-ir,t,k)G^(t,y,k), 
c 

o # = o_ 1 . . ( 31 ) 

Here, the complete connection is given by 

A(— ir, 2°, ft) - A s (-iT,k)5 s (z°) + A a {-iT,k)6 a (z°) 

+ J dv° D{-iT,v Q ,k)I!£{v Q ,z Q ,k) , 
I 

_ _ _ = __l_ + . (32) 

While the first line is already known from the perturba- 
tive case, the second line is a new contribution. It con- 
tains the non-local part of the thermal self-energy. The 
quantities A s , A a and D are straightforward generaliza- 
tions of their perturbative counterparts Ag, Aq and D 
(see eqs. (25,B12) and eqs. (29,B13)). 

In contrast to the perturbative case, eq. (32) is an im- 
plicit equation for the complete connection. For example, 
for the 2PI setting-sun approximation, it has the form, 



+- + 



Equation (32) can be solved iteratively, 

A<°) = A s 5 s (z°) + A a S a {z°) = \— , (33) 

A (*+i) = A s S s (z Q ) + A a S a (z°) + J dv°DI$\ Mk) . 

i 

For example, for the setting-sun approximation, the first 
three steps of this iteration are 



+ — f <£-> 




The first line represents the zeroth step and the first step, 
and the second line shows all diagrams contributing at 
the second step. All diagrams are generated with the 
correct symmetry factors. 

On the other hand, the nonperturbative generalization 
of eq. (29) reads 



G^(-ir,-ir',k) 

= D(-iT,-ir',k) + J dw° J dz° A(-ir, w°, 
c c 

y.G c t C{w\z\k)A T {z\- lT \k) 
= D(—ir, —ir' ,k) 

+ A s (-it, k)G th (0, 0, k)A s (-iT', k) . 



k) 



(34) 



The derivation of this equation is also shown in appendix 
B. Similar to the perturbative case, the formalism estab- 
lished above can be used to relate any Feynman diagram 
formulated on the thermal time path ( "C+X" ) with a set 
of Feynman diagrams formulated on the closed real-time 
path C containing non-local effective vertices represent- 
ing the thermal initial correlations ("C + a"). This is 
accomplished by three steps: 

1. First, the contour integrations over the thermal 
time path C+X associated with internal vertices 
are split into two integrations over C and I. A di- 
agram with V vertices is thus decomposed into 2 V 
contributions. 

2. Next, all internal propagator lines connecting a real 
and an imaginary time are replaced using eq. (31). 
Additionally, the internal propagator lines connect- 
ing two imaginary times are replaced according to 
cq. (34). The parts containing imaginary times are 
encapsulated, which can be visualized by joining 
the complete connections to boxes surrounding the 
imaginary vertices. 
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3. Eventually, the iterative solution of eq. (32) for 
the complete connection is inserted. Each result- 
ing contribution can be identified as a diagram for- 
mulated on the closed real-time path C containing 
non-local effective vertices a n . The latter are con- 
structed explicitly, as appropriate for the underly- 
ing 2PI approximation. 

The first two steps are analogous to the perturbative 
case, with complete propagators and connections instead 
of free ones. The third step is special for the nonpcr- 
turbative case. It results in contributions which contain 
non-local effective vertices a n of arbitrarily high order n. 
These take into account the infinite sequence of thermal 
initial n-point correlations, which are present due to the 
underlying nonperturbative approximation. 

For example, for the nonperturbative setting sun dia- 
gram, step one and two can be written as 



S(x,y) 



-e- 



x 



4 o 



+ • i C^j i • • 

For the second diagram, the third step can be written as 








>- 









X<>- + -K>+<>* + • ■ 

x x-x 



The first diagram on the right-hand side is obtained by 
inserting the zeroth iteration for the four complete con- 
nections. The other diagrams are obtained by inserting 
the first iteration. Note that all diagrams shown above 
are generated with correct symmetry factors. 

Each of the boxes with thin lines represents a non-local 
effective vertex, encoding the correlations of the initial 
state. A thin box that is attached to n propagator lines 
represents a contribution to the thermal initial n-point 
correlation function a„ (xi, . . . , x n ). For example, the 
leading contributions to the thermal initial 4- and 6-point 
correlations are given by 

a 4,0L [Zl,---,Z4) 

= -iX J d 4 v A<°> (v, «i)A(°) {v, z 2 )A ( V ( v , z 3 )A<°) (v, z 4 ) 




X X 
X X 



(35) 



and 

th,2PI, 
a 6.0L \ z 1t ■ ■ ■ i ■ 



= {-iX} 2 J d 4 vj d 4 w A^(v,z 1 )A^(v,z 2 )A^(v,z 3 ) 
x x 

x D(v,w)A^(w,z i )A {0 \w,z 5 )A (Q \w,z 6 ) , 




(36) 



Note that these contributions are nonperturbative ap- 
proximations of the exact thermal initial correlations, 
since they involve the complete thermal propagator. 



V. KADANOFF-BAYM EQUATIONS WITH 
THERMAL INITIAL CORRELATIONS 

In this section, Kadanoff-Baym equations that can de- 
scribe thermal equilibrium on the closed real-time path 
C with finite initial time t init = are derived. This re- 
quires to take into account non-Gaussian correlations of 
the initial (thermal) state of the system. These thermal 
initial correlations have to be determined in accordance 
with the nonperturbative 2PI formalism underlying the 
Kadanoff-Baym equations. Therefore, the techniques de- 
veloped in the previous section are combined with the 
Kadanoff-Baym equations for non-Gaussian initial states 
derived in section II. As before, the main idea is to match 
the Kadanoff-Baym equations describing thermal equi- 
librium within the C+a-formalism on the one hand with 
equivalent evolution equations obtained from the 2PI ef- 
fective action within the C +1- formalism on the other 
hand. 

On the one hand, we use that the Kadanoff-Baym 
equation for a thermal initial state is a special case of 
the Kadanoff-Baym equation for a non-Gaussian initial 
state (see eq. (19)), which has the form 

(^o + k 2 + M 2 th ) G th (x a ,y°, k) = -tS c (x° - y°) 



i J dz° 



+ iU thAa (x ,z°,k) 



G th (z a ,y°,k) . 



(37) 



Here Hg nl (x ,z°,k) and U^ nl (x° , z° , k) denote the 
Gaussian- and non-Gaussian parts of the non-local self- 
energy, respectively. Furthermore, 

U thAa (x°,z°,k) 

= U thtXatF (x°,k)6 s (z°) - l -U thtXatP (x°,k)5 a (z°) 

denotes the contribution from the non-Gaussian initial 
correlations that is supported at the initial time surface 
z° = (see section II). 



12 



On the other hand, the equation of motion of the com- 
plete thermal propagator based on the thermal time con- 
tour ("C+Z") evaluated for x°,y° & C is 



{d 2 x0 + fc 2 + Ml) G th (x°,y , k) = -i5 c+x (x° - y°) 



-/ ! dz°U^(x 
c+x 



(38) 



Of course, the thermal propagator is time-translation in- 
variant and thus only depends on the time difference 
x° — y . The upper notation is chosen for convenience, in 
order to simplify the comparison with the corresponding 
Kadanoff-Baym equations. 

Since the two formulations of thermal equilibrium are 
equivalent, the solutions of the equation of motion (38) 
evaluated for x° 7 y° G C as well as the Kadanoff-Baym 
equation (37) for a thermal initial state have to agree. 
Now we will use this equivalence in order to determine the 
non-Gaussian parts of the thermal self-energy appearing 
on right-hand side of the Kadanoff-Baym equation. 

In order to do so, we start from eq. (38) and split the 
contour integration on the right-hand side into one inte- 
gration over C and one over T. We note that the result- 
ing two contributions are in general not time-translation 
invariant, but their sum certainly is. The integral over 
the imaginary-time contour I can be rewritten using the 
complete connection (32), 

J dz° Ut(x°,z ,k)G th (z ,y°,k) 
x 

<(x°, v°, k) J dz° A(v°,z°, k)G th (z°, y°, k) 
c 

= J dz° J dv°X\t(x { \v { \k)^(v { \z { \k) 

C X 

+ J dw" D(v°, w°, fc)<( W °, z°, fc)) G th (z°, y°,k). 
x 

Inserting this into the equation of motion (38), it takes 
precisely the form of the Kadanoff-Baym equation (37) 
for a non-Gaussian initial state. By comparison, the non- 
Gaussian contributions to the self-energy for the thermal 
initial state can be inferred, 
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(39) 



U^ nl (x°,z°,k)^ J dv a J dw° \lt(x\v\k) 
x I 

x D(v°,w a ,k)U\i(w a ,z a ,k)\ x0 z0eC , 



m th , Xa (x ,z°,k) = I dv" n?» u ,fc) 



,.0 Tinl/J) „.0 



A(°V,z°,fc) 



For the setting-sun approximation, the steps listed above 
leading from the formulation of the Schwinger-Dyson 
equation on the thermal time path ( "C+X" ) to the formu- 
lation on the closed real-time path with thermal initial 
correlations ("C+a") are 
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-0— +-(j>^ + -(j> • j- 



Thus, the Gaussian and non-Gaussian contributions to 
the self-energy in setting-sun approximation for a ther- 
mal initial state are given by 



ng, nl (x°,z°,k) 

iU thAa {x°,z°,k) 




• j : j > 



(40) 



We stress that the propagator determined by the Kada- 
noff-Baym equations for a thermal initial state is time- 
translation invariant, as required. The matching proce- 
dure guarantees that this property is directly inherited 
from the C+I formalism. Nevertheless, the contribu- 
tions to the self-energy derived above are not, and can- 
not be, time-translation invariant individually (compare 
with the discussion in section III and Refs. [45-47]). 

Finally, in order to explicitly obtain the thermal ini- 
tial correlations that are appropriate for a specific 2PI 
approximation, the iterative expansion (33) of the com- 
plete connection has to be inserted. This yields a series 
expansion of the non-Gaussian self-energies, 
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and analogously for n t ^ ; . 

For example, in setting-sun approximation, the ther- 
mal initial correlations obtained from inserting the ze- 
roth, first and second iteration of the complete connec- 
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The zeroth contribution contains the thermal non-local 
effective 4-point vertex (35). The first iteration yields 
three diagrams with thermal effective 6-, 8-, and 10-point 
vertices, and the second iteration yields six contributions 
with thermal effective 8-, 10-, 12- (two diagrams), 14-, 
and 16-point vertices. The smallest and largest are shown 
in the last line of eq. (41). The expansion of n™^; con- 
tains thermal non-local effective vertices of order six and 
higher, 
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(0),nG 



n 



(l),raG 
th.nl 




The zeroth contribution contains the thermal non-local 
effective 6-point vertex (36). The first contribution con- 
tains 15 diagrams with thermal effective vertices of order 
8 to 18. The order of the thermal initial correlations 
appearing up to the fifth contribution in setting-sun ap- 
proximation are shown in table I. 

Thus, vacuum or thermal initial states entail an infi- 
nite hierarchy of initial n-point correlation functions. In 
setting-sun approximation, the non-Gaussian initial cor- 
relation of lowest order is the 4-point correlation given 
in the first line of eq. (41). Its contribution to the Kada- 
noff-Baym equations is 

m th,\ a ( x > z ) = J d4x i23 Gth(x 1 x 1 )G t h(x 1 x 2 ) 
c 

.. r< I \ ■ th,2PI / \ 
x Gth(x, X3)ia it Q L (xi,x 2 ,X3,z) . 

In the limit x°,y° —> 0, only the thermal initial 4-point 
correlation contributes to the right-hand side of the Ka- 
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TABLE I: Thermal initial correlations in 2PI setting-sun ap- 
proximation. The column number is the order n = 4, 6, . . . 
of the thermal initial n-point correlation. The row number 
k = 0, 1, . . . shows which initial correlations contribute to 



^thXa ( u PP er table) and n^'^" G (lower table), respectively. 
Due to the Zn -symmetry, only even correlations are non-zero. 

danoff-Baym equations (37): 

dz° [ng inl + U^ nl + iU thAa } (x°, z°, k)G th (z°, y°, k) 
-> |dz zn^ A J 2 ; ,z ,fc)G 4/l (z ,y ,fc). (42) 



The reason is that all other contributions contain at least 
one internal classical vertex. However, these vertices are 
accompanied by memory integrals. Due to the struc- 
ture of the closed real time-path, memory integrals van- 
ish when all external time arguments approach the initial 
time. 

Altogether, Kadanoff-Baym equations that describe 
thermal equilibrium have been derived on the closed real- 
time path with finite initial time. This requires to take 
into account non-Gaussian n-point correlations of the ini- 
tial state, which enter explicitly on the right-hand side 
of the Kadanoff-Baym equations. The values of these 
initial correlation functions for a thermal state can be 
determined iteratively in accordance with the underlying 
2PI approximation. 
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FIG. 6: Time evolution of the equal-time propagator Gf(£, i,fc) obtained from Kadanoff-Baym equations with thermal initial 
2-point correlation function (initial state (A), red dashed lines) as well as thermal initial 2- and 4-point correlation functions 
(initial state (B), green solid lines), for three momentum modes, respectively. The dotted horizontal lines show the renormalized 
thermal propagator Gth(0,0, k) which serves as initial condition att = 0. 



VI. NUMERICAL RESULTS 



initial 4-point correlation: 



In this section, we support the findings of the pre- 
vious section by numerical solutions of Kadanoff-Baym 
equations. Of course, one cannot implement the com- 
plete infinite hierarchy of initial n-point correlation func- 
tions, which would be required for an exact description 
of thermal equilibrium in the framework of Kadanoff- 
Baym equations. However, thermal equilibrium can be 
approached closer and closer as one includes more and 
more thermal initial correlations. Therefore, we consider 
initial states which are obtained by keeping thermal ini- 
tial ro-point correlation functions for n < n max and set- 
ting all higher initial correlations to zero, 



a n (xi 



5 Xn) 



.Hi 



(Xx, . 





for n < n max , 
for n > n max ■ 



More precisely, we compare the time evolution obtained 
for two ^-symmetric initial states (A) and (B) with 
n-max = 2 and n max = 4, respectively. Both initial states 
have in common that the statistical propagator is ini- 
tialized with the thermal propagator and that all initial 
correlations a n {x\, . . . , x n ) with n > 4 vanish: 



G F (x ,y ,k)\ x o =y0=0 
d x oGF(x°, y , k)\ x o =y o =Q 



Gth(—iT,0,k)\ T ->Q , 
=y"={> — , 

d x od y oG F (x ,y°,k)\ x o =y o =0 = d*Gth(-ir, 0,fc)| T ^ o , 
a n (xx j • • • , x n ) — , for n > 4 

where Gth{—iT, 0, fc) is the complete thermal propagator 
at temperature T = Tmit- The only difference is in the 



a A (xx, ■ 



,x 4 ) = 



{xx 



for (A) , 
, xx) for (B) , 



where a^ 1 (xx, ■ ■ ■ , xx) is chosen as shown in eq. (41). Ac- 
cordingly, (A) is a Gaussian initial state and (B) is a min- 
imal non-Gaussian initial state, and both states "are as 
thermal as possible" for the respective classes of initial 
states. 

In both cases, we employ the 2PI three-loop approxi- 
mation as discussed in Ref. [9] and section II. We note 
that an equivalent set of equations can be obtained from 
the 4PI three-loop approximation [50]. In this case, 
ax determines the initial value of the 4PI 4-point func- 
tion [49]. 

The numerical solutions were obtained on a lattice 
with 32 3 x 2000 2 lattice sites and lattice spacings of 
a s m R = 0.5 and a t m R = 0.025. We use \ R /4\ = 0.75 
for the renormalized coupling constant 4 . The bare mass 
tub and coupling Xb, as well as the thermal propaga- 
tor Gth(—ii~, 0, k) are determined by a separate numerics 
program as described in Refs. [34, 35] . 



4 We checked for our lattice settings that this value for the cou- 
pling is far below the 'Landau pole' exhibited by the 2PI 3-loop 
approximation [35]. Apart from that, we also checked that the 
damping rates for the unequal-time propagators increase like A 2 , 
which indicates that the above value for the coupling lies in the 
domain of validity of the 2PI 3-loop approximation. 
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FIG. 7: Time evolution of the effective temperature and effective chemical potential obtained from Kadanoff-Baym equations 
with thermal initial 2-point correlation function (initial state (A), red dashed lines) as well as thermal initial 2- and 4-point 
correlation functions (initial state (B), green solid lines). The shaded areas illustrate qualitatively the deviation of the effective 
particle number density n(t, k) from the Bose-Einstein distribution function. They are obtained from the asymptotic standard 
error of the fit (via least-square method) magnified by a factor 10, for better visibility. Nevertheless, the errors become invisibly 
small at times tmR S> 10. 



A. Time-evolution of the equal-time propagator 

Both initial states (A) and (B) are nonequilibrium ini- 
tial states, since the thermal initial correlations higher 
than Umax are neglected. However, among all possible 
Gaussian initial states, the state (A) has the minimal 
offset from thermal equilibrium. Similarly, the state (B) 
has the minimal offset from equilibrium among all possi- 
ble initial states that are parameterized by initial 2- and 
4-point correlation functions. In the following, we will 
compare this minimal offset for both initial states (A) 
and (B). In exact thermal equilibrium, the equal-time 
propagator Gp{t,t,k) is constant. Therefore, the time- 
dependence of the equal-time propagator is a measure for 
the deviation from thermal equilibrium. 

Figure 6 shows the time-evolution of the equal-time 
propagator obtained from the numerical solution of the 
two sets of Kadanoff-Baym equations. We find that the 
deviation from thermal equilibrium is considerable for 
the Kadanoff-Baym equations with Gaussian initial state 
(A). In contrast to this, for the Kadanoff-Baym equa- 
tions with the non-Gaussian initial state (B) the equal- 
time propagator remains very close to the thermal prop- 
agator at all times. We have checked that this qualita- 
tive behaviour stays the same when varying the coupling 
strength or the lattice spacings. This shows that the ther- 
mal initial 4-point correlation already yields a reasonable 
approximation of the complete thermal initial state. 

This observation can be understood by analyzing the 
role of the thermal initial correlations (40) in the Kada- 
noff-Baym equations (20). In fact, there are two distinct 
reasons why the contribution of the thermal initial n- 



point correlations are suppressed for n > 4. 

The first reason is that the effective loss of memory 
is stronger the larger n. This can be seen as follows: 
The thermal initial n-point correlations enter the Kada- 
noff-Baym equations as effective n-point vertices. These 
are connected with classical vertices by n lines. Each of 
these lines yields a propagator G(t, 0, fe) for which one 
of the time arguments is evaluated at the initial time 
(see eq. (9)). However, such unequal-time propagators 
are damped exponentially with respect to the time t [5]. 
Therefore the contribution of the initial n-point correla- 
tions is also damped exponentially for t — > oo, and the 
damping is the stronger the larger n. 

The second reason is that the contributions to the Ka- 
danoff-Baym equations (20) which arise from thermal ini- 
tial n-point correlations with n > 4, vanish in the limit 
t — > as has been shown in section V. 

Thus, the influence of thermal initial n-point corre- 
lations with n > 4 on the solutions of Kadanoff-Baym 
equations is suppressed compared to the initial 4-point 
correlation for early times (t — ► 0) as well as for late 
times (t — * oo). 

B. Offset between initial and final temperature 

Figure 7 shows the time evolution of the effective tem- 
perature T(t) and the effective chemical potential 5 fi(t) . 



5 Note that the effective chemical potential is introduced here as a 
fit-parameter in order to characterize the equilibration process. 
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They are obtained by fitting a Bosc-Einstcin distribution 
function to the effective particle number density n(t, k) 
for all times t. The effective particle number density can 
directly be extracted from solutions of Kadanoff-Baym 
equations [9]. At the initial time, we have T(0) = T ini t 
and /i(0) = 0. For t — > oo, the solutions approach ther- 
mal equilibrium, i.e. T(t) — ► Tfi na i and — ► within 
the numerical accuracy. The time-evolution can be di- 
vided into the three phases of (i) correlation build-up for 
tmR < 1, (ii) kinetic equilibration for 1 < tm R < 10 
and (hi) chemical equilibration for tm R > 10. The nu- 
merical solutions clearly exhibit a separation between the 
time-scales of kinetic and chemical equilibration [9, 22]. 

We see in figure 7 that the deviation from thermal 
equilibrium is much smaller for the non-Gaussian initial 
state (B) as compared to the Gaussian initial state (A). 
This result supports the observation from the previous 
subsection. 

In addition, figure 7 reveals a qualitative difference be- 
tween both sets of equations. For the Gaussian initial 
condition (A), there is an offset between the initial tem- 
perature Ti n i t an d the final temperature Tf ina i- However, 
for the non-Gaussian initial condition (B), the initial and 
final value of the temperature agree within the numerical 
accuracy. 

This behaviour can be explained by the composition of 
the total energy of the system. It can be split into kinetic 
and correlation energy, 



E total — Ekin(t) + E corr (t) , 



(43) 



where 
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+ 



1 fd 3 k 



2 J (2tt) 3 

A rd?q 

4 .fe) 



d x od y o + k 2 + m 2 



/' d 3 q 
^-^G F (t,t,q) G F (x°,y ,k)\ x o =y o = t , 



E r 



-(*) - - 



d 3 k 



4 J (2tt) £ 
1, 



(n Xa . F (t,k)G F (o,t,k) 



+ -n Xa , p (t,k)G p (o,t,k) 



z 
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The total energy is conserved by the numerical solutions 
up to numerical errors (< 1%). Using the thermal ini- 
tial state, we can also derive an expression for the total 
energy in thermal equilibrium at temperature T, 

E e \T) = EH(T)+Etl rr {T) , 



Since particle number is not conserved for A<E> 4 -theory, the effec- 
tive chemical potential has to vanish in thermal equilibrium. 



^here E%f n (T) = E km (t = 0)\ G=Gth and 
1 f d 3 k 

Ecorr(T) = ~- J-^^ flth, Xa, F(t, k)G t h(0,t,k)\t=0 ■ 

For all possible initial states, the final temperature can 
then be determined by the requirement 

Etotal = E cq (Tfi na i) ■ 

Both initial states considered here feature thermal n- 
point correlations at temperature T n it for n < n max . 
Therefore, we have 

Etotal — E eq (T init ) — AE( n >n max )(Tinit) j 

where AE^ n >n max )(T) denotes the contribution to the 
thermal energy that comes from n-point correlations with 
n > Umax- For initial state (A), we find 

AE (n>2) (T) = EZAT) . 

Since the thermal correlation energy has a non-zero value, 
we find that 

Tfinai ^ Tnit for initial state (A). 

This is in accordance with the results of Ref. [28] for 
the non-relativistic case. In general, one would expect 
that this is also true for initial state (B). However, using 
cq. (42), we obtain 



E eq 

corr 



(T) = * / ^ 



1234 G t h(x,Xi)G t h(x,X 2 ) 



x G t h(x,x 3 )ial (xi,X2,X3,Xi)Gth(xi,x)\ x =o ■ 

Thus, only the initial 4-point correlation contributes 
to the thermal correlation energy. Therefore we have 
Ai?( n> 4) (T) = 0, which means that the initial and final 
temperatures have to agree, 



- final 



— Tj/t 



for initial state (B). 



This property is quite remarkable. It means that the to- 
tal energies for the initial state (B) and for the complete 
thermal initial state are identical. Thus, the thermal ini- 
tial n-point correlations o?^ with n > 4 do not contribute 
to the total energy of the initial state. Instead, the ther- 
mal initial 4-point correlation af 1 already captures the 
complete thermal correlation energy at the initial time 
in setting-sun approximation. 

Moreover, this result can be generalized to arbitrary 
2PI approximations. In general, the correlation energy 
at the initial time is given by 



E corr (t — 0) 



d 3 k 

4 7(2^)3 



U Xa , F (t,k)G F (0,t,k)\t=o • 



The most general structure of the non-Gaussian self- 
energy Ii\ a (x,y) is shown in figure 8. The Kernels 
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FIG. 8: Contribution Tl\ a (x,y) to the self-energy Tl(x,y) 
where the left line is connected to a classical vertex, and the 
right line to an effective non-local vertex. 



A^ and A^ G may in general contain classical as well 
as non-local effective vertices. In setting-sun approxima- 
tion (13), they are given by A2° — on and A^ G = 0. 
In the limit t = x — > 0, all contributions containing 
classical vertices vanish due to the memory integrals ac- 
companying these vertices. Thus, at the initial time, only 
those parts of A™ G and A^ G do contribute to the correla- 
tion energy that can be rewritten in the form of an initial 
4- and 3-point vertex, respectively. 

This property is characteristic for the ^-interaction. 
If the Lagrangian would contain a (non-renormalizable) 
$ 6 -interaction, then also initial 5- and 6-point correla- 
tions would explicitly contribute to the energy density of 
the thermal initial state. 



VII. CONCLUSIONS AND OUTLOOK 

In this work, we derive Kadanoff-Baym equations for 
non-Gaussian correlated initial states describing vacuum 
and thermal equilibrium. We provide suitable techniques 
to establish approximations to the exact thermal ini- 
tial correlation functions that match the approximation 
scheme underlying the Kadanoff-Baym equations. These 
techniques are applicable for arbitrary truncations of the 
2PI effective action. Examples are given for the 2PI 
three-loop approximation. 

Finally, we discuss numerical solutions of Kadanoff- 
Baym equations for a real scalar <E> 4 quantum field the- 
ory, which take the thermal initial 4-point correlation as 
the leading non-Gaussian correction into account. These 
solutions are compared to solutions obtained for Gaus- 
sian initial states. For the latter, the initial state has no 
correlation energy by definition. Therefore, even if one 
initializes the two-point function with the thermal prop- 
agator for a certain temperature, the system equilibrates 
at a different temperature [28]. We show numerically 
and analytically that this feature of the Gaussian ini- 
tial state is remedied completely already by taking the 
thermal initial 4-point correlation into account. The rea- 
son is that higher correlations of the initial state do not 
contribute to the total energy at the initial time. Thus, 
we find that including an initial 4-point correlation func- 
tion yields a significantly improved approximation to the 
complete thermal initial state as compared to Gaussian 
initial states. 

The techniques developed in this work provide a frame- 
work for investigating the renormalization of Kadanoff- 
Baym equations. It is known that the 2PI effective ac- 
tion can be renormalized in thermal equilibrium [34, 35] . 



Accordingly, it is possible to derive renormalized Kada- 
noff-Baym equations for the thermal initial state using 
the techniques introduced above. These provide a well- 
defined expansion point for nonequilibrium initial states. 
For example, it is possible to parameterize the initial n- 
point correlation functions in the form a n = a„ + Aa„ . 
In the limit Aa n — > 0, the thermal initial state is recov- 
ered and thus the Kadanoff-Baym equations are formally 
finite. Then it is possible to investigate which conditions 
the deviations Aa„ have to fulfill such that no new di- 
vergences are introduced. This is left to future work. 
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APPENDIX A: PERTURB ATIVE CALCULATION 
OF THERMAL DENSITY MATRIX ELEMENT 

In order to describe the equilibrium limit within non- 
equilibrium quantum field theory, it is important to cal- 
culate the thermal density matrix element {<p+ \pth \ <P—} 
evaluated with respect to two eigenstates $(Q,x)\<p±) = 
(p±(x)\cp±) of the field operator $(f, x) at time t — ti n u = 
0. In the following, the perturbative expansion of the 
thermal density matrix element is discussed based on 
Rcf. [44]. Therefore, the action (1) is formulated on 
the imaginary time path, and split into the free part 
£o[<fo] = f x d 4 x ((<9(/>o) 2 — m 2 <t>Q) /2 and the interaction 
part Sint[4>] = — Jx^ 4x I\ < l } ( x ) 4 - For the free thermal 
density matrix Po = exp (—(3H ) containing the free 
Hamiltonian Hq, which is quadratic in the field, the ma- 
trix element can be calculated by a path integral that is 
analogous to eq. (21). The result is [44] 



(<P+ \Po\ <P-) = A/" exp iS [(j)o} 



where Mo is a normalization factor, which is independent 
of tp± , and (f>o (x) is the solution of the free equation of 
motion dSo/S(f> = (— □ — m 2 )(j) = 0onl subject to the 
boundary conditions 

0o (0, x) = (p-(x) and c/) (-i[3,x) = ip + (x) . 

The solution is uniquely determined, and, in spatial mo- 
mentum space, given by 

sinh(w fc r) sinh(w fe (/3 - r)) 
(fio(-iT, k) = -—- -^(f+ik) H —7 ^ <f-( k ) . 



sinh(wfc/3) 
where oJi = m, 2 + k 2 . 



sinh(wfe/3) 
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The full thermal initial correlations can be obtained by 
perturbing the full Hamiltonian H around Hq, 



(<P+ \Pth\<P-) =A/"exp i(S [<p ] +F mt [<p }) 



where TV is a normalization factor and i<So[0o] is the free 
contribution. iF int {<j>o] is the sum of all connected Feyn- 
man diagrams with vertices given by the derivatives of 
SintW\ evaluated for <f>(x) = 0, 



int 



8<tA 



= -i\5x{xi-x 2 )5j:(xi-x?,)5z{xi-x i ) = 



According to the Feynman rules given in section III the 
empty circle reminds us that the corresponding integra- 
tion runs over the imaginary-time contour I. The bound- 
ary conditions of the path integral (21) are formally taken 
into account by the "field expectation value" 

M~ir,k) = (g), 

along the imaginary contour T, as well as the propagator 

= (Al) 

- ^sinhK/?)! 8 ^^^^^'^^-^ 
+ sinh(w fe r') sinh(cj fe (/3 - r))9(r - r'fj , 

which is the Greens function for solutions of the free equa- 
tion of motion that vanish at the boundaries r = 0, (5. Dq 
is denoted by the dotted line. To first order in A, iF int [4>o] 
is given by 



D (-it, -it', k) 
1 



intl<P0\\ O ( X ) 



+ 



4! 



d 4 x [3Dl(x,x)+6(f)l(x)D Q (x,x)+(t>%(x) 



The field-independent diagrams, like the first one above, 
can be absorbed into the normalization TV. The per- 
turbative expansions of the thermal initial correlations 
a% are obtained by the n-th functional derivative of 
F[ip + , ip_] = Sq[4>o} + Fi n t{(j)o\ with respect to the field, 



in 



th,e! 



' ( x l > • • • j x n) 



5iF[if + ,Lp- 



SifeAx!) . ..5tp €n (x n ) 



ip=0 



to which all diagrams with n insertions of 4>o contribute. 
Here, the decomposition from eq. (9) has been used. 
Thus, the initial correlations obtained in this way are 
indeed supported only at the initial time, as required. 
Formally, the functional derivative corresponds to replac- 
ing the field insertions by (distinguishable) external lines 
in the diagrammatic expansion of iF int [<po] according to 

<p (-iT,k) i ► A (-ir, x°,k) 



•I- 



where Ao(— it, x°, k) is defined in eq. (24). For example, 
the leading contribution to the thermal initial four-point 
correlation function obtained from the fourth derivative 
of iF int [<fi } is given by (see eq. (27)) 



X x 

:o: 

X" X 



ial h 0L {xi,X2,x 3 ,x 4: ) = V (A2) 

= -i\ J d 4 v A (v,x 1 )A (v,x 2 )A (v,x 3 )A (v,x i ) . 



Switching again to momentum space, an explicit expres- 
sion for the leading contribution to the perturbative ther- 
mal initial four-point correlation function is obtained, 



• t/l,ei,€2|€3i«4 



i a 



4,0L 



where 



p 

(k 1 ,k 2 ,k 3 ,k 4 ,) = -A J d,T Aq 1 (-fr, fei 



xAq 2 {-it, k 2 )A% (-it, fc 8 )A§* {-it, fc 4 ) , 



A+(-iT,k) = 

Ao(-ir,fc) = 



sinh(wfcr) 
sinh(w fe /3) 
sinh(q> fc (/3 - t)) 
sinh(wfc^) 



The integral over the imaginary time can be performed 
analytically. In the zero-temperature limit (/3 — ► oo), one 
obtains 



la 



4,0L 



(fci,fc 2 ,fe 3 ,fc4) 



(A3) 



-A 



for ei = e 2 = e 3 = e 4 = ± , 



__ / ^fc x + UJ k 2 +^fc 3 +^fc 4 

\ else . 

Altogether, a diagrammatic expansion of the matrix ele- 
ment of the thermal density matrix in terms of perturba- 
tive Feynman diagrams has been developed as suggested 
in Ref. [44]. This allows to explicitly calculate thermal 
correlation functions order by order in the quartic cou- 
pling constant. The lowest-order perturbative result (A2) 
may be compared to the nonperturbative 2PI result (35). 



APPENDIX B: DERIVATION OF THE 
CONNECTION IN THE 2PI CASE 

In this appendix, we generalize eq. (23) to the 2PI for- 
malism. Before we consider the complete thermal prop- 
agator, we find it helpful to take an intermediate step 
by first considering a mixed thermal propagator which 
is identical to the complete thermal propagator on the 
imaginary branch of the thermal time contour, and which 
obeys the free equation of motion on the real branches. 
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1. Mixed thermal propagator 

It is straightforward to define projectors on the parts 
C and X of the thermal time contour, 



and 



lx(z°) 



lc(z°) = 



if x° e c 

1 if 



1 if x° G C 



v if x° eX ' 

Of course, they fulfill the relation 

li(x°) + l c (x°) = 1 for all x° G C+X . 

The mixed thermal propagator is defined by the following 
equation of motion, 

G m,th( x >V) = i( n x + m 2 )S c+x (x - y) 

-l x (x°)l x (y°)IL th (x,y), (Bl) 

where x° 7 y° G C+X. Here U t h(x,y) is the complete 
thermal self-energy. The mixed propagator can be de- 
composed into statistical and spectral components, 

G m ,th(x,y) = G m , F (x,y)~^sign c+x (x a -y°)G m ^ p (x,y) . 

The equation of motion for the mixed propagator can 
equivalently be written as 

(p x + m 2 ) G m>t h(x, y) = -iS c +x(x - y) 

- il x {x Q ) j ' d 4 z U m ,th(x, z)G mtth (z,y) . 
x 

Each of the two time arguments of the propagator can ei- 
ther be real or imaginary, which yields four combinations 

G %m G %th' G * C th> G ^th- Thc mixcd propagator eval- 
uated with two imaginary time arguments is identical to 
the complete thermal 2PI propagator, 

G nuth(x, V) = G ¥(x, y) for x°, y° G X, (B2) 

whereas the mixed propagators evaluated with one or two 
real time arguments, G^ tt (a;,y) and G^ th (x,y), fulfill 
the equation of motion of the free propagator, 

{n x +m 2 )G% th (x,y) = 0, (B3) 
(u x +m 2 )G% F (x,y) = {u x + m 2 )G% p (x,y) = 0. 

At the initial time x° = y° = 0, thc propagators on all 
branches of the thermal time path agree. Using eq. (B2), 
one obtains 

G ^Ah( x >y)\x°= y °=o = G th{x,y)\ x o= y o =0 , 



for Vi G {C,X}. Thus, the initial value of the mixed 
propagator at x° = y = is given by the complete 
thermal propagator. 

For the mixed propagator with one imaginary and one 
real time, the equation of motion, transformed to spatial 
momentum space, reads 

(-a T 2 + fc 2 + m 2 )G^Hr,y°,fc) 



- J dr' U^(-ir,-ir',k)Gl e th (-iT',y ,k) . 



Next, a Fourier transformation with respect to the imag- 
inary time is performed, using in particular 



J dr cxp(-iuj n T)d 2 Gf^ th (-iT,y a ,k) 
o 

= -ul G ^th(Un,y°,k) 

+ disc(iu; n Gl c th + d T Gl c th )(y°, k) , 

where ui n = 2ir/3n is a Matsubara frequency. It is impor- 
tant to take the contribution from boundary terms into 
account, 

disc^G^+^G^J^fc) 

= [(iu; n GZ% + d T GZ%)(-iT,y ,k)} T T 2o • 

Thus, the Fourier transformed equation for thc mixed 
propagator reads 

(c 2 +fe 2 + m 2 )Gf ra c f , l K, 2 ; ,fc) (B4) 
= -nf*(u n ,k)Gl c th (u; n ,y ,k) 

+ disc( lw „ G^ th + d T G^ th )(y a ,k) . 

The boundary terms have to fulfill the equation of motion 

(a 2 + fe 2 + m 2 ) disc(Gl c th )(y°, k) = , 
(5 2 o+fe 2 + m 2 )disc(9 r G^)( 2 /°,fe) = 0, 

which follows using G^ c f/l (w„, y°, fe) = G^ th (y° ,u n ,k) 
and the equation of motion (B3) for G^ th . Furthermore, 
the initial conditions at y° — are fixed by the period- 
icity relation of the thermal propagator as well as the 
equal-time commutation relations, 

disc(G* c tt )(0,fc) = 0, 

d y o disc(G^ c fft )(0,fc) = i , 

disc(d r G* c t/l )(0,fc) = 1, 

d y adi S c(d T G^ th )(0,k) = 0. 

The statistical and spectral components G^ c F (0,y°,fc) 
and G%f (0,y°,k) of the mixed propagator are two lin- 
early independent solutions of the free equation of mo- 
tion. Since it is a second order differential equation, any 
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solution can be expressed as a linear combination, espe- 
cially 

disc(Gl c th )(y°,k) = -iG% p (0,y°,k), 

G£>(o,y°,fc) 



disc(d T Gl e th )(y°,k) = 



Gth(P,0,k) 



Inserting this result into eq. (B4) and using the Fourier- 
transformed Schwinger-Dyson equation (30) for the com- 
plete thermal propagator yields 



Gl c th (u n ,y°,k) 



Gf n c F (o,y°,k) 



2. Complete thermal propagator 

Using the equation of motion (Bl) of the mixed prop- 
agator, the self-consistent equation of motion (30) of the 
complete propagator can be rewritten as 

G^{x,y) = G-] th (x,y) - [l - l x (x°)l x ( y ")} n th (x,y) , 

for x°,y° £ C+X. By convolving this equation with Gth 
from the left and with G m .th from the right, the integral 
form of the Schwinger-Dyson equation is obtained: 



(0,0,*). 
(iu n Gf?(uj n ,k)) G% p (0,y°,k). 



Finally, the upper relation can be rewritten in the form 
Gl e th (u; n ,y ,k) = J dz° A m K, z°, k)G% th (z°, y°, k) , 

(B5) 

where a mixed connection has been introduced, 
A m (w n , z°, k) = A^K, k)6 s (z°) + A^K, k)6 a (z°) 



G t h(x,y) = G m ,th(x,y) + J du j d v G th {x,u) 

c+x c+x 

x [1 - l x (u°)l x (v )] ILth(u,v)G m ,thC",y) ■ (B7) 

Evaluating it for x° C and y G X, and performing a 
Fourier transformation with respect to the relative spa- 
tial coordinate x — y as well as the imaginary time y° 
gives 

G%(x°,u n ,k) = G e * th (x°,w n ,k) 



{ cliOAk) ) ^ + y 2iuJ nGfh(^n,k)J S a (z°) + j du oj dv o Gth{x ^ u ^ k) 



1— . (B6) 

Furthermore, the transposed connection is defined as 
A^ n (z°,uj n ,k) = A m (uj n , z°,k). Eq. (B5) for the mixed 
propagator is the generalization of eq. (23) for the free 
propagator. Thus, the mixed propagator evaluated with 
one real and one imaginary time can be written as the 
convolution of the mixed connection, which involves the 
complete 2PI propagator, and the real-real mixed prop- 
agator, which obeys the free equation of motion. 



c+x c 



xU th (u ,v°,k)G^; th (v ,uj n ,k) 
-i J du° G%(x°,u°,k)Il th (u ,u; n ,k)Gl x th (u> n ,k) 



Next, G^ th (x°,u> n ,k) and G^ th (v°, u> n , k) are replaced 
using eq. (B5) with interchanged arguments. Further- 
more, it is used that Gf^ th (uj ni k) = Gf^(u) n ,k) (sec 
eq. (B2)): 
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G%(x°,u n ,k) = J dz° G% h {x ,z\k)A T m {z°,uj n ,k) 
c 

+ Jdz° J du° J dv° G th (x ,u a ,k)n th (u°,v^k)Gi%(v^z^k)Al(z n ,u n ,k) 

C C+I c 

-i I du° G%(x°,u°,k)Il th (u ,w n ,k)Gf?(co n ,k) 

G%(x°,z°,k)-J du°J dv° (G^(x ,u ,k)U th (u ,v ,k)Gl c th (v ,z°,k) 



c 

dz° 



-i J du a G%(x°,u°,k)Tl th (u ,w n ,k)Gf?(uj n ,k) 
c 

J dz a G%(x°,z°,k) |A^(z°, w „,fc) -zn f ,(z°, w „,fc)Gi :r K,fc) 
-J du° J dv° U th (z ,v a ,k)Gl c th (v Q ,u ,k)Al(u a ,uj n ,k) 



A^(*V n ,fc) 



In the second step, the Schwinger-Dyson equation (B7) additional term, which is the convolution of the thermal 

evaluated for x°,z° G C was used again. In the third self-energy, evaluated with one real time and one Matsub- 

step the complete real-real propagator was factored out axa, frequency, with, the propagator D^uj mi u} ni k). This 

by interchanging the integration variables u° z a in the propagator is given by 
second and third term. The last line can be simplified by ^ 

Fourier transforming with respect to the imaginary time D(uj n , Lu m , k) = Gf/f (u> n , k) 
v°, and performing the integral over C using eq. (B6): 

.0 „,o i.\r<TC /,o „,o . U T/„o 



du v J dv» Ti th {z»y, k)G^ th {v", u", fc)A^( u u , w „, k) 

C X 

= -iTj2^th(z°, u>i, fe)G^(o; ; ,0, k)A s m (w n , k) 



dw° dz" A(w„, w°, k)Gth(w u , z°, fc) A 1 u m , fc) 



c c 

5n ' m -Gff(w n , k) - G *h^k)G^{Lo m ,k) (B9) 



T 



Gtfc(0,0,fc) 



= -*r^n tft (z°,^,fc)Gi :r (^,fc) 



Gfh"(Vn, fc) 

Gth(0,0,fc) 



In the last line 

/dz° X(w„, W °, fc)G tt (u; , z°, fe)n a (z°, w m , fc) = 



Finally, a decomposition of the complete thermal 2PI wag ^ where ^ £ {A ,U th }. The propagator D has 
propagator evaluated with one real time and one Mat- ^ e properties 



subara frequency is obtained, 

G%(x°,u; n ,k) = J dz° G%(x ,z°,k)A T (z°,w n ,k) , 
c 

where the complete connection was introduced, 

A T (z°,u Jn ,k)^Al(z ,u Jn ,k) (B8) 
-iT^2 Kth{z ,u) m ,k)D(u m ,oJn,k) , 

m 

with A(u> n ,z°,k) = A T (z°, ui n , fc). Compared to the 
mixed connection, the complete connection contains an 



D((j n ,oJ m ,k) = D(u} m ,w n ,k) , 
T^D{u n ,w m ,k) - 0. 

m 

From the last property it can be inferred that only the 
non-local part of the thermal self-energy U t h(z a , u m , fc) = 
^■th + TIth( z °> u m, fc) contributes in eq. (B8), since the 
local part is independent of the Matsubara frequency. 

By applying an inverse Fourier transformation with re- 
spect to imaginary time, using in particular 

D(-ir, -it', k)=T 2 J2 «m, fc) , 
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the complete thermal 2PI propagator with one imaginary 
and one real time can be decomposed as 

Cff(x°,-iT,k) = Jdz»G c t £{x*,z*,k)A T {z*,^T,k), 



Eqs. (B10,B11,B12) constitute the nonperturbative gen- 
eralizations of eqs. (23,24). 

The nonperturbative generalization of eq. (29) is ob- 
tained from eq. (B9), 



and 



Gf h c (-ir,y°,k) = Jdz° A{-iT,z\k)G c t £{z\y\k) 



o • 



f 



where the complete connection is given by 
A(-ir, z°, k) = A m (-ir, z°, k) 

+ J dv° D(-ir,v ,k)U^(v ,z°,k) 
i 

= A s (-ir, k)S s (z°) + A a (-ir, k)5 a (z°) 
+ J dv° D(-iT,v°,k)U^(v ,z°,k) 



+ — (n nl 



(B10) 



(Bll) 
and 

A T (z°, — it, fe) = A(-ir, z°, k) = — | . 

The coefficients A s,a (— ir, k) are given by 

( ' ' G th (0,0,k) ' 
A a {-iT,k) = 2d T Gf*(-iT,0,k) . (B12) 



Gf h x (-ir, -ir',k) 



D{-iT,-ir',k)+ J dw° J dz° A(-ir, w°,k) 



c c 

„0 ,0 maT/J „• / 



xG t/l ( W °,z",fc)A J (z°,-ir',k) 
D(—ir, —ir', fe) 

+ A s Ht, k)Gth(0, 0, fc)A s (-ir', fe) . (B13) 



Note that only the parts of the connections containing 
A s contribute to the double integral in the second and 
third line, whereas the parts involving A a and n™^ vanish 
due to a cancellation of the contributions from the two 
branches of the closed real-time path. Using the Feyn- 
man rules from above, the upper equation can also be 
written as 



-o + o- 



I 1 



In Summary, there are two differences compared to 
the perturbative case: (i) the free thermal propaga- 
tor Go,th{— iT, 0, fe) enters the free connection, whereas 
the complete thermal propagator Gth{— it, 0, fe) enters 
the complete connection, and (ii) the free connection 
Ao(— it, z°, fe) is only supported at the initial time z a = 
0±, whereas the complete connection A(— it, z°,k) fea- 
tures an additional term containing the non-local part of 
the complete thermal self-energy. 
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